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We use kinetic Monte Carlo simulations to understand growth- and etching-induced step bunching 
of 6H-SiC{0001} vicinal surfaces oriented towards < 1100 > and < 1120 >. By taking account of the 
different rates of surface diffusion on three inequivalent terraces, we reproduce the experimentally 
observed tendency for single bilayer height steps to bunch into half unit cell height steps. By taking 
account of the different mobilities of steps with different structures, we reproduce the experimentally 
observed tendency for adjacent pairs of half unit cell height steps to bunch into full unit cell height 
steps. A prediction of our simulations is that growth- induced and etching- induced step bunching 
lead to different surface terminations for the exposed terraces when full unit cell height steps are 
present. 

PACS numbers: 81.10.Aj, 81.15.Kk, 81.65.Cf 



I. INTRODUCTION 

Silicon carbide (SiC) is a very promising material for 
microelectronic applications because of its superior elec- 
tronic properties, high thermal and chemical stability, 
high-power and high-frequency capability, and high tol- 
erance to radiation damagei*^ SiC is also an attractive 
candidate as a substrate for the heteroepitaxial growth of 
other materials. A particulary exciting example (which 
motivated the present study) is the growth of epitax- 
ial graphene by thermal decomposition of the basal sur- 
faces of single crystal 4H and 6H SiC- Nevertheless, SiC 
will not reach its anticipated potential until a variety of 
problems are solved, not least being the need to con- 
trollably grow device-quality single crystal material on a 
large scale ii^i^ 

One approach to the growth problem is "step- 
controlled" epitaxy, where new layers grown onto sur- 
faces vicinal to the hexagonal basal planes inherit the 
stacking order of the substrate through the step-flow 
mode of growth.^ Unfortunately, step-flow growth on 
vicinal surfaces does not always proceed by the uni- 
form motion of a train of evenly spaced steps. Instead, 
growth-induced step bunching often occurs, as it in- 
variably does when vicinal surfaces are etched by expo- 
sure to hot hydrogen gas. Suggestions for the origin of 
SiC step bunching include impurity adsorption, differ- 
ences in surface energetics for different bilayers of a-SiC 
poly types ^ii^iiiii^ differences in intrinsic step velocities 
and step configurations, ^ i -'^^i and other differences in 
growth kinetics. ^^^^^ However, no systematic exploration 
of any particular mechanism and comparison of the re- 
sults with all available data seems to have been performed 
until now. 

This paper reports the results of kinetic Monte Carlo 
(KMC) simulations designed to identify the kinetic path- 
ways that promote growth- and etching-induced step 
bunching of vicinal 6H-SiC surfaces. We focus on sur- 
faces vicinal to (0001) (Si-terminated) and (OOOT) (C- 



terminated) with steps running perpendicular to the 
< iTOO > and < 1120 > directions. These par- 
ticular starting surfaces were chosen to make con- 
tact with experimental observations made using atomic 
force microscopy (AFM), low-energy electron diffrac- 
tion (LEED) analysis, high-resolution transmission elec- 
tron microscopy (HRTEM), and scanning tunnelling mi- 
croscopy (STM)j^ii^iii^^^'^^^^^^^^ Our main conclusion is 
that the experimental results for the Si-face are quite ex- 
plicable using a lattice model that recognizes that there 
are three inequivalent terraces for surface diffusion and 
two inequivalent steps with different mobilities. The C- 
face data are similarly explicable (or at least rationaliz- 
able) if the terrace diffusion rates and step mobilities are 
less different on this face than on the Si-face. 



II. PREVIOUS WORK 

A. Experimental Observations 

This section reviews the experimental observations of 
etching- and growth-induced bunching for surfaces vici- 
nal to the basal planes of 6H-SiC. The designation 6H 
refers to the bulk unit cell indicated in Figure [1] where 
six bilayers of silicon and carbons atoms are arranged in 
the particular stacking sequence shown. If it were ex- 
posed to vacuum, the top layer of silicon atoms would 
be a typical (0001) Si-terminated surface. The bottom 
layer of carbon atoms similarly exposed would be a typ- 
ical (0001) C-terminated surface. Along the c-axis of 
6H-SiC crystal, the bilayers are arranged in two groups 
of three bilayers each. Each bilayer is exactly the same 
apart from rigid lateral shifts within one group of bilayers 
and rotation of the lattice by 60° around the c-axis from 
one group of bilayers (A, B, C) to the other {A* ,C* ,B*). 
The hexagonal arrangement of atoms shown in Figured 
is a (projected) view along [0001] of one of the bilayers 
from the 74*(7*5*-group seen in edge-view in Figure [H 
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FIG. 1: Arrangement of Si and C atoms in 6H-SiC crystal. A 
fragment of (1120) atomic plain is shown schematically. 
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TABLE I: Experimentally observed etching-induced 
step bunching of (0001) (Si-terminated) and (OOOT) (C- 
terminated) surfaces of 6H-SiC for various miscut angles and 
orientations. Entries in quotation marks are inferred by the 
present authors. See text for discussion. 
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FIG. 2: Top view of a single bilayer is shown schematically. 
Li and Li are the units of length in < 1100 > and < 1120 > 
directions, correspondingly. 



erly be regarded as "nanofacets" .^^ For the same surface 
miscut 12° along < 1120 >, the authors find no aver- 
age step orientation due to large-scale step meandering. 
From the corresponding image of the C-face, we infer 
(hence the quotation marks in the table) similar, but less 
pronounced, step meandering. Quite generally, the data 
summarized Table I demonstrate that the tendency for 
etching-induced bunching is greater on the Si-face than 
on the C-face and greater for steps oriented perpendicu- 
lar to < 1100 > than for steps oriented perpendicular to 
< 1120 >. 



There have been many experimental studies of 
growtbi^^OaOi^^^ and etchingi^O^Oa^^ of 6H SiC 

surfaces cut slightly vicinal to its basal planes. The de- 
tails vary (temperature, doping, and Si/C ratio during 
growth), but most papers focus on the Si-terminated face 
with miscuts oriented along the < 1120 > direction. Only 
a few report results for the C-terminated surface or for 
vicinal miscuts oriented along the < 1100 > direction. 
For our purposes the systematic H-etching experiments 
conducted by Feenstra and co-workers are particularly 
valuable.^^ Table I summarizes their AFM and HRTEM 
observations for nominally on-axis samples and samples 
with intentional miscuts of ~ 3° and 12°. The numbers 
I and 1 in the table refer to observations of more-or-less 
ordered arrays of bunches where 3 Si-C bilayer steps have 
bunched into a single step with the height of \ unit cell 
and where two such bunches have further bunched into 
one unit cell height steps. Etching of (0001) surfaces 
miscut 12° along < 1100 > produces increased bunch- 
ing into 4-5 unit cell height steps which may more prop- 



There is no single data set for growth-induced step 
bunching comparable to the etching-induced results sum- 
marized in Table I. Nevertheless, a survey of the liter- 
ature reveals trends very similar to the etching data. 
Thus, for surfaces miscut by 3.5° toward < 1120 >, 
Kimoto et al? find that growth on Si-terminated sur- 
faces produces | unit cell height bunches while growth 
on comparable C-terminated surfaces is twice as likely 
to remain completely unbunched (only SiC bilayer steps 
appear) as to bunch into \ unit cell height steps. Sim- 
ilarly, data obtained for step-flow growth on vicinal Si- 
terminated surfaces oriented toward < 1100 > exhibit 
6 bilayer bunches (full unit cell), compared to only 3 bi- 
layer bunches (half unit cell) observed for similar surfaces 
miscut along < 1120 Therefore, as during etching, 
the tendency for growth-induced bunching is greater on 
the Si-face than on the C-face and greater for steps ori- 
ented perpendicular to < 1100 > than for steps oriented 
perpendicular to < 1120 >. For both growth and etch- 
ing, step bunching is always more pronounced on surfaces 
with higher miscut angles. 
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B. KMC Simulations 

Three groups have used kinetic Monte Carlo simula- 
tions to study step flow growth on vicinal SiC{0001} sur- 
faces. Heuell used a one-dimensional model that did not 
distinguish carbon atoms from silicon atoms Two in- 
equivalent types of steps were considered and the proba- 
bilities for a diffusing adatom to attach to each step from 
the terraces below and above were treated as independent 
parameters. A parameter set was found were an initial 
train of height-one steps bunched into a train of height- 
six steps. However, trains of | unit cell height steps are 
commonly observed in experiments, which suggests that 
a train of single bilayer steps bunches first into a train of 
I unit cell height steps, which then bunch into full unit 
cell height stepsji^i^U Heuell's model does not produce this 
behavior. 

Stout developed a very elaborate KMC simulation that 
took account of the SiC crystal structure, the transport, 
adsorption, and surface diffusion of physisorbed precur- 
sors, the dissociative chemisorption, surface diffusion, 
and desorption of dissociated species, and the attach- 
ment/detachment of adatoms to/from step edges The 
energy barrier for a particular atom to make an activated 
Monte Carlo move was taken to be proportional to the 
product of the coordination numbers of the initial state 
site and the final state site. For growth onto surfaces 
vicinal to 6H SiC{0001}, Stout's simulations bunched an 
initial train of single bilayer steps into a train where two 
nearby single bilayer steps accompany a 4 bilayer height 
step. No bunching into 3 bilayer height steps or 6 bilayer 
height steps was observed. 

Finally, Camarda and co-workers used a full lattice 
KMC model including defect sites to study step flow 
growth onto surfaces vicinal to 4H SiC(OOOl).^^ These 
authors did not treat silicon and carbon atoms as diffus- 
ing species; the smallest growth unit considered was a 
Si-C dimer. Bunching was observed, but the step heights 
were not reported. 



III. SIMULATION MODEL 

We have developed a three-dimensional KMC simu- 
lation model based on the crystal structure of SiC. In 
this paper, the model is used to study step-flow etching 
and step-flow growth of surfaces vicinal to 6H-SiC{0001}. 
Later work will address island nucleation and multilayer 
roughness on singular surfaces and thermal decompo- 
sition of stepped and flat surfaces to produce epitax- 
ial graphene. The starting vicinal surface studied was 
usually a uniform train of 36 single bilayer steps with a 
miscut angle of ~ 15° (25°) for miscut oriented towards 
< 1120 >(< iTOO >) . Otherwise (see Figure [2j), miscuts 
oriented towards [1100] were treated using "helicoidal" 
boundary conditions (HBC)^^ along the [1010] direction 
and periodic boundary conditions (PBC) along the [1120] 
direction. The typical system size was 216Li x 40^2. 



Miscuts oriented towards [1210] were treated using HBC 
along the [1120] direction and PBC along the [lOTO] di- 
rection. The typical system size was 72Li x 2I6L2. 

The standard KMC method^^ identifies a set of ele- 
mentary "moves" and catalogs their relative rates. For 
our simulations, thermal desorption directly into the gas 
phase was not allowed, but all atoms (except fully coor- 
dinated bulk atoms) were permitted to move to empty 
nearest-neighbor or next-nearest-neighbor surface sites 
(with equal probability) at a rate R = Rq ex.p{—E/kBT), 
where Rq = 10^^ /sec, /c^ is Boltzmann's constant, and 
T is the substrate temperature ( T :^ 1000 K for most 
of our simulations). The activation energy E depends 
on the atom type and its local coordination through a 
bond-counting rule that includes only the four possible 
nearest neighbors. All the simulations we report used 
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with £;si-c = 0.75 eV, £;si-si = 0.35 eV, and Ec-c = 
0.65 eV. The absolute values of these parameters are not 
crucial because they only represent effective pair-bond 
energies. What matters is their relative ordering, which 
reflects (i) the stability of the SiC crystal and (ii) the 
much greater strength of the C-C bond compared to the 
Si-Si bond. 

We come now to the crucial feature that distinguishes 
our simulations from others. For one species type (C or 
Si), the surface jump rate computed using Eq. (1) is ex- 
actly the same when the atom sits on any of the six 6H 
(OOOl)-type terraces (called A, B, C, A*, C*, and 5* 
in Figure [1]) exposed by a vicinal surface with only sin- 
gle bilayer steps. However, the beyond nearest-neighbor 
interactions that energetically distinguish the many dif- 
ferent polytypes of SiC from one other imply that the 
energy barriers to surface migration cannot be exactly 
the same for all six terraces. The relevant surface dif- 
fusion barriers have not been reliably computed or mea- 
sured. However, the first-principles, surface total energy 
calculations of Righi et al.^^ show a clear energetic pref- 
erence for SiC (0001) surfaces to continue their subsurface 
stacking order. This conclusion agrees with the observa- 
tion that 3C polytype islands nucleate on 4H- and 6H- 
SiC(OOOl) substrates and with observed stable surface 
terminations for these exposed faces. ^^^^^'^^ Accordin gly, 
we use the surface energy ordering computed in Ref. |28| 
to scale the energy barriers for terrace surface diffusion. 
The scaling factors used were 1.0 for the atoms sitting 
on A and A* terraces, 1.15 for the atoms sitting on C 
and 5* terraces, and 1.3 for the atoms sitting on B and 
C* terraces. This procedure is consistent with our use of 
the binding energies in Eq. (1) to estimate the bare en- 
ergy barrier. Our choice of the scaling factors, or, more 
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FIG. 3: Growth-induced step bunching on vicinal 6H- 
SiC(OOOl) surface with miscut towards [1100]. 



FIG. 4: Etching-induced step bunching on vicinal 6H- 
SiC(OOOl) surface with miscut towards [1100]. 



exactly, their ordering, is in agreement with the rela- 
tive stability of the inequivalent terraces implied by the 
three different step velocities observed during the step 
flow growth of graphene by the decomposition of vicinal 
6H-SiC(0001).3i 

Etching was simulated very simply. Every time an 
atom was selected to move along the surface, it was in- 
stead removed entirely from the simulation with a proba- 
bility e, where < e < 1. The etching results were quite 
insensitive to the exact value of e. To study growth, we 



ignored precursor effects and deposited silicon and car- 
bon with equal probability at randomly chosen empty 
surface sites of the SiC lattice. The average deposi- 
tion rate used, F = 100/s, corresponds to a typical SiC 
growth rate of ~ 1 fim/h. 



IV. RESULTS 

Our model produces very similar results for growth 
and etching of the Si-terminated face and C-terminated 
face of 6H SiC. This differs from the experimental obser- 
vations summarized in Section 11. For that reason, this 
section reports results only for the Si-terminated face. 
The C-terminated face is discussed in Section V(D). 

A. Si-face with vicinal miscut towards < 1100 > 

Figure [3] shows a sequence of simulated morphologies 
during step flow growth onto a vicinal 6H-SiC(0001) sur- 
face with the miscut towards the [1100] direction. The 
step bunches that eventually form have the height of one 
6H unit cell (6 bilayers). Note, however, that ^ unit 
cell height steps (3 bilayers) form first. The faster grow- 
ing bunch catches up with the slower growing 
bunch (ABC) to form the final full unit cell height bunch. 
The step-edges of the final bunches are mostly straight 
and aligned along [1120] direction, which is the energet- 
ically most stable configuration. The overhanging step 
risers produced by the simulation is an artifact of the 
nearest-neighbor approximation used in (1). A less steep 
(and smoother) step occurs when next-nearest-neighbors 
are included, but at the cost of a much slower simula- 
tion. This change introduces no qualitative effects on 
the bunching, so we used only the simpler model in this 
paper. 

For comparison with Figure [3l Figure [4] shows a se- 
quence of simulated morphologies during the step flow 
etching of a vicinal 6H-SiC(0001) surface with the mis- 
cut towards [1100] direction. The etching morphology we 
obtain is similar to the growth morphology except that 
the (A*C*5*) bunch retracts faster than (ABC), and 
ends up at the bottom of the full unit cell height step. 
During growth, the {A*C*B*) bunch winds up on top of 
the (ABC) bunch. It is interesting to note that many ex- 
periments precede epitaxial growth of 6H-SiC with a gas 
etching step to smoothen the surface. To study this case, 
we performed a growth simulation beginning with the 
surface shown in the last panel of Figure [H As shown 
in Figure [5l the starting full unit cell height step with 
(ABC) on top of {A*C*B*) flips during the growth to a 
fuh unit ceh height step with (A*C*5*) on top of (ABC). 
In other words, the starting surface does not matter and 
the final panel of Figure [3] is the same as the final panel 
of Figure [H 
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FIG. 5: Growth- induced step bunching on vicinal 6H- 
SiC(OOOl) surface with miscut towards [1100]. Initial surface 
configuration corresponds to etching-induced train of one unit 
cell height steps. 



B. Si-face with vicinal miscut towards < 1120 > 

Figure [6] shows a sequence of simulated morphologies 
during step flow growth onto a vicinal 6II-SiC(0001) sur- 
face with the miscut towards the [1210] direction. We 
observe the formation of ^ unit cell height steps with 
zigzag shaped step edges. On average, the steps are 



aligned along the [1010] direction (perpendicular to the 
miscut direction). However, the straight segments of the 
step edges are aligned along the close-packed < 1120 > 
directions. Etching produces very similar results in the 
sense that ^ unit cell height steps form. This is shown 
in Figure [3 



[0001] 




FIG. 6: Growth-induced step bunching on vicinal 6H- 
SiC(OOOl) surface with miscut towards [1120]. 



V. DISCUSSION 

The results of our KMC simulations are in good qual- 
itative agreement with experiments. The formation of 
full unit cell height steps (six bilayers) is inherent in SiC 
step-flow growth/etching on vicinal 6H-SiC(0001) sur- 
faces with the miscut towards < 1100 >. On the other 
hand, ^ unit cell height steps form during simulations of 
growth/etching when the miscut is along < 1120 >. We 
understand all these features in terms of the six differ- 
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FIG. 7: Etching-induced step bunching on vicinal 6H- 
SiC(OOOl) surface with miscut towards [1120]. 



ent steps that appear on any surface composed of only 
bilayer height steps (top panel of Figure [H or Figure [7j) . 
There are six different velocities because there are three 
inequivalent terraces where surface diffusion events occur 
and two inequivalent step edges, where attachment, de- 
tachment, and interlayer transport events occur. As de- 
scribed in Section III, the surface diffusion rate is fastest 
for adatoms sitting on A and A* terraces, slower for 
adatoms sitting on C and B* terraces, and slowest for 
adatoms sitting on B and C* terraces. Figure [8] shows the 
two types of steps. ^^^^^ For the Sn step, a next-to-next 
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FIG. 8: Side view of Sn and Sd steps. See text for discussion. 
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nearest neighbor jump is required for an atom to attach 
to the step from its upper bounding terrace. Moreover, 
after this jump occurs, the attached atom is only singly 
bonded to the step, and thus easily detached. For the 
Sd step, only a next-nearest neighbor jump is required 
for an atom to attach to the step from its upper bound- 
ing terrace. After this jump occurs, the attached atom is 
doubly bonded to the step, and thus less likely to detach. 

Figure [9] and Figure [11] show scenarios for growth and 
etching, respectively. The fastest steps during growth are 
bounded from above by the most stable terraces (A and 
A*) . The slowest steps during growth are bounded from 
above by the least stable terraces (C and B*). For the 
same reason of stability, fast growing steps are the slowest 
etching and vice versa. Accordingly, the terraces exposed 
at the time when | unit cell height steps are present on 
the surface, in both cases correspond to A and A* (see 
Figure [1]), which is in agreement with experiments^ 



A. Step-Flow Growth 

The top panel of Figure [9] shows a vicinal surface com- 
posed of single bilayer steps. The arrows on the steps re- 
flect their relative velocities due to the rates of surface dif- 
fusion on the three inequivalent terraces mentioned just 
above. As a result, the single bilayer steps bunch into 
the ^ unit cell height steps shown in the second panel of 
Figure [9l The subsequent bunching of these steps into 
the full unit cell height step shown in the final panel of 
Figure [9] occurs because one half cell height step is Sn 
type and the other is Sd type. The latter moves faster 
than the former because the next-to-next nearest neigh- 
bor jump required for attachment to the Sn step from 
above is not allowed in our model. In reality, we presume 
there is simply a higher barrier for this process to occur 
at a Sn step than at a Sd step. In other words, the 
Ehrlich-Schwoebel barriers associated with these steps 
are different. ^^^^^ 

The preceding discussion applies directly for vicinal 
miscuts toward the < 1100 > direction. However, we 
have stated that full unit cell height steps do not form 
in our SiC growth simulations for miscuts toward the 
< 1120 > direction. This occurs because the step edges 
in this case have a natural zigzag shape, consisting of 
alternating straight segments, corresponding to Sn^ or 
Sd steps (see Figure [TQ|). Each step edge has equal por- 
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FIG. 9: A cartoon of step- flow growth on 6H SiC(OOOl): (a) 
Different length arrows indicate the different growth velocities 
of the steps which terminate the three inequivalent bilayer 
terraces; (b) Two types of | unit cell height steps {Sn and Sd) 
differ by the number of dangling bonds for the outermost step 
edge atoms. The presence of an Ehrlich-Schwoebel energy 
barrier to downward inter layer diffusion at Sn steps explains 
the difference in growth speed between Sn and Sd ^ unit cell 
height step bunches; (c) Sd steps wind up on top Sn steps 
to form a single, full unit cell height step. 



tions of Sn and Sd steps, and, as a result, all ^ unit cell 
height steps propagate with the same speed. Of course, if 
the miscut is not exactly toward the < 1120 > direction, 
adjacent | unit cell height steps differ in their relative 
population of Sn and Sd step edges. This may trigger 
the formation of full unit cell height steps. 
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FIG. 10: (a) A step edge of the so called "open" step (top 
view), aligned along < 1100 >, is shown schematically, (b) 
Growth and etching at such steps typically results in develop- 
ment of triangular protrusions. The alternating straight seg- 
ments of this protrusions {Sn and Sd step edges) are aligned 
along the < 1120 > directions. 
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B. Step Flow Etching 



Etching is often regarded as the inverse of growth^i^A^S 
Therefore, mimicking our discussion of growth, the top 
panel of Figure [11] shows a vicinal surface composed of 
single bilayer steps. The arrows on the steps reflect their 
etch velocities due to the different step detachment rates 
associated with the three inequivalent terraces discussed 
earlier. As a result, the steps bunch to form the ^ unit 
cell height steps shown in the second panel of Figure [Til 
The subsequent bunching of these steps into the full unit 
cell height step shown in the final panel of Figure [9] oc- 
curs because the ^Ar-type step, which has triply bonded 
outermost atoms, etches more slowly than the ^'^-type 
step, which has only doubly bonded outermost atoms. 

In agreement with experiments, we observe the devel- 
opment of "triangular" protrusions (see Figure [4|) which 
form as a result of etching at Sd step bunches. As Fig- 
ure [12] shows, the straight segments of these protrusions 
are aligned at angles of 30° with respect to the direc- 
tion of miscut ([1100]) and thus correspond to energeti- 
cally stable Sn step bunches. As discussed in detail in 
Ref. |33|, the outermost atoms of these protrusions, which 
have only two bonds with nearest neighbors, constitute 
another source of instability. As a result, preferential 
etching of these protrusions leads to the formation of full 
unit cell height steps. 

In our simulations of SiC etching for miscut towards 
< 1120 >, we do not observe the formation of unit cell 
height steps. This agrees with the most recent exper- 
imental results We explain this in terms of the 
previously discussed zigzag structure of the steps which 
occur on this surface. Each ^ unit cell height step with 
a zigzag shape has equal portions of faster and slower 
etched straight segments, which correspond to Sd and 
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FIG. 11: A cartoon of step-flow etching of 6H SiC(OOOl): (a) 
Different length arrows indicate the different etch velocities 
of the steps which terminate the three inequivalent bilayer 
terraces; (b) Two types of | unit cell height steps (Sn and 
Sd) have different edge velocities because they differ in the 
number of dangling bonds for the outermost step edge atoms; 
(c) Sn steps wind up on top Sd steps to form a single, full 
unit cell height step. 



Sn step bunches. For this reason the etching rates of all 
I unit cell height step are identical (on average) and unit 
cell height steps do not form. 

Notwithstanding the foregoing, some etching experi- 
ments on surfaces with miscut towards < 1120 > do see 
the formation of full unit cell height steps. ^^^^^^^^ A pos- 
sible explanation for these conflicting observations is the 
previously mentioned possibility of deviations of the mis- 
cut from exactly < 1120 > with its attendant steps with 
faster and slower etching segments. When the popula- 
tions of these segments is not equal, the | unit cell height 
steps with more fast-etching segments catch up to the 
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FIG. 12: (a) Sd step edge (top view) is shown schematically, 
(b) Sd steps are less stable then Sn steps, because the out- 
ermost atoms of the former have two nearest-neighbor bonds, 
while the outermost atoms of the latter have three such bonds. 
Growth and etching is faster at Sd steps and typically results 
in development of triangular protrusions. The straight seg- 
ments of this protrusions are aligned along the < 1120 > 
directions and correspond to energetically stable Sn steps. 



steps with fewer fast-etching segments. The result is 
a train of full unit cell height steps. 



C. Surface Terminations 

Only one type of terrace is exposed to the vacuum 
after all full unit cell steps have formed on a [1100] mis- 
cut surface. However, our KMC simulations predict that 
the exposed terrace is not the same after growth as it 
is after etching. After growth, the A* C*B* -hunch (the 
outermost atoms of bilayer-steps have two bonds with 
nearest neighbors and two dangling bonds) is on top of 
the ABC-bunch (the outermost atoms of bilayer-steps 
have three bonds with nearest neighbors and one dan- 
gling bond). This is called the surface termination 
in the literature. After etching, the sequence of bilay- 
ers at the surface is opposite: ...B*C* A*CBA. This is 
called the surface termination. Cross section TEM 
experiments could be used to test this prediction. 



for the C-face compared to the Si-face. ^^^^ On the other 
hand, the C-face data are quite explicable if the terrace 
diffusion rates and step mobilities differ from their values 
on the Si-face. 

The terrace diffusion scaling factors used to take ac- 
count of the three inequivalent terraces of SiC(OOOl) were 
chosen based on the Si- face calculations of Righi et. al?^^ 
These authors did not perform similar calculations for 
the C-face and it is possible that the results are differ- 
ent. One possibility is that the scale factor ordering is the 
same as for the Si-face, but that their magnitudes are less 
different. Another possibility is that the ordering of the 
scale factors differs on the C-face. Experimental support 
for this comes from the different surface termination ob- 
served for the two polar faces of 6H-SiC{0001}. The so- 
called (2 X 2)^ reconstruction, which stabilizes hexagonal 
stacking at the surface ( Si stable surface termination) is 
sometimes observed on the C-face. This implies a dif- 
ferent ordering for the three inequivalent terrace scaling 
factors than for the Si-face. 

It seems quite likely that the Ehrlich-Schwoebel bar- 
riers to interlayer diffusion (which strongly influence the 
bunching of half unit cell height steps into full unit cell 
height steps in our model) are different on the C-face 
and the Si-face. The magnitude of these barriers is inti- 
mately connected to the structure of the steps and there 
is theoretical evidence that the step structure indeed de- 
pends on the polarity of the surface*^ If the difference 
between the corresponding barriers is less pronounced for 
the C-face, we would find the experimentally observed de- 
lay in the formation of full unit cell height steps during 
the growth of 6H-SiC on the vicinal 6II-SiC(0001) sur- 
face. Moreover, strong barriers to interlayer diffusion at 
both Sn and Sd steps suppress interlayer mass transport 
and stabilize the persistence of single bilayer steps. This 
would explain the observed experimentally preference of 
bilayer height steps on C-terminated surfaces to remain 
completely unbunched during growth.^ 

Finally, it is possible that step bunching is less pro- 
nounced on the C-face because the elastic-driven repul- 
sive interaction between the steps of a vicinal surface^^ is 
more pronounced on the C-face compared to the Si-face. 
The repulsion depends on the step stiffness, which in 
turn depends on the step structure, which are doubtless 
different for the two faces. 



D. The C-terminated surface 

The simulation results we have presented so far de- 
scribe the evolution of surface morphology during the 
epitaxial growth/etching of 6H-SiC on the vicinal 611- 
SiC(OOOl) surface (Si-terminated face). Our simulation 
results for the C-terminated face (6II-SiC(0001)), using 
the same model parameters, are qualitatively very simi- 
lar. This disagrees with the step bunching behavior ob- 
served in experiments which is typically less pronounced 



VI. CONCLUSIONS 

We have used lattice KMC simulations to study the 
formation of step bunches during growth and etching 
of 6II-SiC(0001) vicinal surfaces. For both situations, 
the simulations show that single bilayer steps bunch into 
half unit cell steps (3 bilayers each), which subsequently 
bunch into full unit cell steps. This is consistent with ex- 
perimental observations for both the Si-terminated face 
and the C-terminated face except that we obtain greater 
bunching for the C-face than seen in experiment. The 
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main driving force for bunching into half unit ceh height 
steps is that surface diffusion is not equahy fast on ah bi- 
layer terraces. The main driving force for the subsequent 
bunching into fuh height unit cehs is the existence of two 
different local atomic step structures, which leads to two 
different step mobilities. A prediction of the model which 
invites an experimental test is that growth-induced and 
etching-induced step bunching lead to different surface 
terminations for the exposed terraces when full unit cell 
steps are present. 
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